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The ti asymmetry measured at the Tevatron continues to disagree with Standard Model 
predictions at the 3 sigma level. We update the status of the phenomenological light axigluon 
model in explaining the asymmetry data, taking into account constraints from the charge 
asymmetry at the LHC and the ti cross section at both Tevatron and LHC. We find that 
an axigluon with a mass between 100 and 400 GeV provides an excellent fit to the data. 
Recent searches by ATLAS and CMS for pair production of heavy resonances which decay 
to dijets rule out axigluons with large branching fractions to dijets. However axigluons 
which predominantly decay to multijets via intermediate resonances are still a possibility. 
We outline four distinct scenarios which cover the most important decay topologies and 
discuss how one might exclude or discover axigluons as multijet resonances at the LHC. 
MadGraph implementations for each of the scenarios are provided. 
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I. INTRODUCTION 

The forward-backward asymmetry A FB as measured in top quark pair production at the Teva- 
tron continues to disagree with QCD predictions Q3431. The asymmetry has been seen in events 
where only one of the top quarks decays leptonically and in events where both top and anti-top 
decay leptonically. Combining the single-lepton DO analysis [5J of 5.4 fb _1 of data with CDF's 
single-lepton analysis (6) of 8.7 fb _1 and CDF's dilepton analysis [7] of 5.3 fb _1 we obtain 1 
A FB = 0.185 ± 0.037 for the "unfolded" (parton level) asymmetry. This is 3.2 a larger than the 
NLO QCD and electroweak prediction of A FB = 0.066 obtained in |6l using POWHEG. Focusing 
on the dependence of the asymmetry on the invariant mass of the top pairs or their distribution in 
rapidity [0[TOl one obtains 3 a deviations without the need to combine experiments. Altogether 
it appears very unlikely that the it asymmetry at the Tevatron is due to statistical fluctuations. 
Standard Model (SM) predictions from different NLO event generators vary by amounts which 
are much smaller than the experimental errors, and NNLO as well as nonperturbative corrections 
are expected to be smaller yet. However, since the asymmetry first arises at NLO a reliable under- 
standing of the theory errors requires a NNLO calculation which is currently in progress. 

Axigluons |fTTT - l20ll with a mass in the range from 100 to 400 GeV can explain the asymmetry. 
Such light axigluons appear to be the "last man standing" 2 of the large number of models (see 
e.g. Il2"5t - l2"7ll for an overview of the models and an extensive list of the original references). Other 
models have difficulty obtaining a large asymmetry at the Tevatron while remaining consistent 
with constraints from the total ti cross sections at the Tevatron and LHC, from the invariant mass 
spectrum of the ti cross sections, from bounds on single top and same-sign top production at the 
LHC, from the dijet cross section at the Tevatron, and from precision low-energy measurements 
such as atomic parity violation ||28l . 

In this paper we first review the current status of ti data in the context of a phenomenological 
light axigluon model (in Section [ll]). We find that the model fits all it data very well. This en- 
courages us to take the model more seriously and discuss constraints from other experiments. The 

1 We used the method of [8 1 for combining the results and ignored correlations in systematic errors by adding statistic 
and systematic errors in quadrature. 

2 Recently, the authors of ||2T1|221 asserted that a complex Z' [23 1 coupling to up- and top-quarks can also explain 
the ti asymmetry and is consistent with all other constraints. The high energy tail of ti production predicted by 
such Z' models differs significantly from SM predictions, but it is consistent with current LHC measurements E4l 
which have large statistical uncertainties. 
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most significant constraints come from pair production of axigluons at the LHC as we show in 



Section III Assuming that axigluons predominantly decay to quark-antiquark pairs one obtains a 
4-jet final state with two dijet resonances of equal mass. This final state has been searched for by 
both ATLAS fl29l ED and CMS [EH and axigluons in the entire mass range from 100-400 GeV 
are ruled out. The beauty of this search is that it is independent of the strength and flavor structure 
of the axigluon coupling to fermions. It only assumes that axigluons are produced via QCD from 
initial state gluons (as required by QCD) and that they decay predominantly to quarks of any of the 
5 light flavors (the searches did not require or exclude b-quarks). Axigluons of mass > 2m t could 
also predominantly decay to top quark pairs if the coupling to top quarks is enhanced. However 
the axigluon pair production cross section is so large (10s of pb for a 400 GeV axigluon at the 7 
TeV LHC) that the resulting 4 top quark final state would have been seen, for example in CMS' 
same sign top search Il3~2ll3~3~ll . 

Thus to save the light axigluon explanation of the tt asymmetry we must postulate that axi- 
gluons do not predominantly decay to dijets but instead decay to a new final state which has not 
been looked for or is very difficult to distinguish from backgrounds. Axigluons carry color, there- 
fore the decay necessarily involves jets. In fact, it should involve only jets because any other 
particles in the final state (leptons, photons, missing Et) would make the signal too easy to detect. 
And since decays to dijets are ruled out by the above-mentioned searches we are led to consider 
models in which the axigluon decays to three or more jets. This can only dominate over the dijet 
decays if axigluons first undergo a two-body decay to intermediate resonances which then further 
decay to the final jets. The existence of such intermediate resonances with large couplings to axi- 
gluons is not far-fetched in axigluon models with small couplings to quarks. In fact, to obtain 
these small couplings one must introduce additional vector-like quarks for the SM quarks to mix 
with, and in the limit where the SM quarks couple weakly to the axigluon the new quarks couple 
strongly. Thus light axigluon models naturally contain additional particles which can catalyze the 
axigluon decays to multi-jet final states. The details of the spectrum of the new particles are model 
dependent and it is interesting to look for the most important possibilities. 



In Section IV we define four consistent axigluon models which cover the most important axi- 
gluon decay topologies. We outline their multijet signatures and discuss current bounds. We find 
that there are several possibilities for a large multi-jet signal to be observed in the near future. To 
facilitate the study of our models we provide ll34l MadGraph ll35l implementations created with 
FeynRules (35|. 
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II. PHENOMENOLOGICAL LIGHT AXIGLUON MODEL AND tt DATA 

The relevant parameters of our phenomenological axigluon model for the ti asymmetry are the 
axigluon mass m A and width Ta and the products of its couplings to the light quarks (gy,g A ) 
and to the top quark (g v , g A ). In the center of mass frame and at tree level, the partonic qq — > tt 
differential cross section including the interference between gluon and axigluon is 



dcr(qq — >■ tt) (5 
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where /3 = ^/l — Am 2 /s. The asymmetric term proportional to c = /3 cos 6 in the interference of 
the axigluon with the gluon (second line) requires non- vanishing axial couplings g A ,g A . Vectorial 
couplings predominantly contribute to the symmetric part of the cross section and are strongly con- 
strained by the good agreement of the measured cross section with SM predictions. We therefore 
choose purely axial couplings to both light and heavy quarks. 3 The axigluon width is important 
in determining the observability of the axigluon as a resonance in dijets at the Tevatron or UA1 
and UA2. For ti production the width is relevant only when m A > 2m t so that top quarks can 
be produced on resonance. Since we are interested in axigluon masses both below and above the 
ti threshold we choose a large width T A = 0.1m A for all of our reference points. The remaining 
parameters which determine the ti cross section and asymmetry at the Tevatron (and LHC) are 
then the product of the axial couplings to light and heavy quarks a A = g q A g A /An and the axigluon 
mass m A . 

In Figure 1 we show the predicted New Physics (NP) contribution to the ti asymmetry Ap§ 
at the Tevatron as a function of a A for 4 representative axigluon masses m A = 100, 200, 300, 400 
GeV. The asymmetry is linear in a A as long as Ap§ ^ 15%, and large asymmetries can be ob- 
tained for moderate values of the couplings. A good fit (1 sigma) to the Tevatron data requires a 
NP contribution to the asymmetry between 8 and 16%. The asymmetry was calculated using Mad- 



3 This restriction, gy — g f v — 0, can be relaxed. Small vectorial couplings do not significantly change the mul- 
tijet phenomenology which is the main focus of this paper. Large vectorial couplings are disfavored by the good 
agreement of tt cross section shape with SM predictions as well as precision fits l28ll37ll . 
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FIG. 1: New Physics contribution to the parton-level forward-backward asymmetry Afb at the Tevatron as 
a function of aA for four different axigluon masses. Also shown are the one and two sigma bands obtained 
by combining the CDF and DO measurements Afb = 0.185 ± 0.037 and subtracting a SM contribution of 
AfAj = 0.066. 

Graph/MadEvent 5, as were all further calculations in the remainder of the paper unless otherwise 
stated. 

At the LHC it production exhibits a small forward-central asymmetry, usually called "charge 
asymmetry" A c . In the axigluon model Ap^ at the Tevatron and A^ p at the LHC are linear 
in aA and therefore the predictions for the two asymmetries are also linearly related for small 
enough aA- Figure 2 shows this correlation for four different axigluon masses. Also shown are 
lines indicating the boundaries of the la preferred values for the asymmetries as measured at the 
Tevatron and LHC. For the LHC number we combined the CMS single-lepton measurement (with 
4.9 fb" 1 ) OH and the combined single-lepton (1.04 fb" 1 ) plus dilepton (4.71 ftr 1 ) result from 
ATLAS [39] . The colored areas correspond to the 68% and 95% preferred regions from our fit to 
the combined LHC and Tevatron measurements. The plot shows that axigluons with any mass in 
the range considered are consistent with asymmetry data at the la level for appropriately chosen 
a A . 

For further comparison with other data we choose four reference axigluon models with masses 
100, 200, 300, 400 GeV and choose axigluon couplings aA (see Table [I]) which produce 10% 
asymmetry from NP at the Tevatron. These models can now be tested against other it data from the 
Tevatron and LHC. Columns 4 and 5 of the Table show the NP contributions to the it cross section 
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FIG. 2: The LHC charge asymmetry and the Tevatron forward-backward asymmetry predicted in four axi- 
gluon models. The vertical and horizontal lines correspond to the the 1 a boundaries of the experimentally 
preferred asymmetries for the Tevatron and LHC (see text). The shaded regions are preferred at 68% (green) 
and 95% (yellow) confidence level in a fit to the combined Tevatron and LHC measurements. 

at the Tevatron and the LHC respectively. In all cases, the new contributions to the cross sections 
are much smaller than the experimental uncertainties of the cross section measurements of about 
±0.5 pb for the Tevatron [@0l|4l] and ±5 pb for the LHC Il42l444ll . Note that since our axigluon is 
light, the cross section enhancement is almost universal over the full range of ti invariant masses 
so that the shape of the cross section da /dM t i does not give interesting constraints on the model. 
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TABLE I: Axigluon coupling strength and NP contributions to top physics in four axigluon models. Each 
model produces a NP contribution to the ti asymmetry of 10% at the Tevatron. The Table gives the NP 
contribution to the LHC charge asymmetry, Tevatron ti cross section and LHC ti cross section. 

In order to increase its discriminating power, CMS also measured the charge asymmetry binned 
by the rapidity of the ti center of mass \y t j\ . In Figure [3] we show the prediction for this differential 
asymmetry due to a 200 GeV axigluon added to the NLO Standard Model contribution compared 
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FIG. 3: Differential charge asymmetry at the LHC as a function of the ti center of mass rapidity (data with 
1 a errors taken from [45 ]). The green dashed line corresponds to the QCD prediction at NLO, and the red 
dotted line is our prediction including the contribution of the 200 GeV axigluon from Table [I] 

with the experimental result [|45l . There is a substantial increase in the predicted asymmetry for 
high \y t i\ which can be understood from the higher percentile of qq initial states in this kinematical 
region. The error bars are still too large to conclude anything definite, but given the large difference 
between the central value of the experimental result and the prediction, the charge asymmetry at 
high y t i will become a very interesting discriminator ll46T - l48l for our model in the future. 

Finally, two additional independent ti physics observables which are sensitive to NP are two 
distinct FB asymmetries of the leptons produced in top decays. The FB asymmetry of leptons 
from events with large ti invariant mass is sensitive to the chirality of the produced top quarks 
ll49ll . whereas the FB asymmetry of leptons from ti pairs produced near threshold is sensitive to 
NP with chiral couplings to the initial quarks in the colliding protons ||501I5T|. The light axigluon 
model is in good agreement with the Tevatron data for both ll6l[T0l. 

III. CONSTRAINTS FROM AXIGLUON PAIR PRODUCTION AT THE LHC 

In this Section we consider recent constraints from ATLAS and CMS which looked for pair- 
production of heavy resonances with subsequent decays to pairs of dijets. As we will see, these 
searches are very powerful in the case of the axigluon because the axigluon pair production cross 
section is enormous. It is enhanced by color and spin factors and is significantly larger than — for 
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example — the pair production cross section of squarks, quarks, or even gluinos of the same mass. 
The cross section is largely model independent because it is dominated by gluon-gluon scattering 
gg — > AA (see Fig. [5]) which is uniquely fixed by gauge invariance and unitarity. 4 
The Lagrangian describing the relevant couplings of the axigluon to gluons is 

C = -hr {D^A U - D V A„) 2 + m 2 A tr (A^) + % X g s tr (G^[A„, A v ]) , (2) 

where D^A U = d^A u — ig s [G^, A v \, A^ is the axigluon field, g s the strong coupling constant, 
the gluon and G^ u is the gluon field strength. 

Notice that beside the usual kinetic term for a colored vector field and the axigluon mass, there 
is another renormalizable operator coupling gluons and axigluons given by the third term. Gauge 
invariance under ordinary color allows an arbitrary value for its coefficient x, however unitarity 
of axigluon pair production amplitudes requires \ — 1 Il52l . To understand this, note that the 
calculation of the amplitude for production of massive vector bosons involves terms which grow 
with the scattering energy. In a consistent unitary theory these terms cancel due to the underlying 
spontaneously broken gauge invariance. This cancellation requires the coefficient x — 1- The 
presence of the x-term can also be seen very easily in a weakly coupled UV completion of the 
axigluon model with the gauge symmetry breaking pattern SU (3)i x SU(3) 2 — > SU (3) co ; or . Here 
the x-term simply arises from rewriting the gauge boson kinetic terms of the mass eigenstates. 
More generally, deviations from x = 1 are consistent with unitarity if they arise at the loop level 
or in the presence of additional massive colored particles contributing to the gg — > A A process. 
We discuss the constraints imposed by unitarity on consistent axigluon models in Appendix|Cj For 
the remainder of the paper we will use x = 1 for our plots. 

Interestingly, while the x-term suppresses the unphysical growth of scattering amplitudes in 
the UV, it enhances the amplitude near threshold where most of the cross section lies. Thus the 
total axigluon pair production cross section at the LHC is considerably enhanced by the x-term 
as can be seen in Figure |4} This cross section is very large. For comparison, the cross section for 
axigluon pairs is only about factor of 5 below the total QCD dijet cross section (with a cut on jet 
p T equal to m^). 

Also shown in Figure [4] are bounds on the axigluon pair production cross section from three 

4 In addition to initial state gluons, there is also a contribution from quark anti-quark collisions (see Figure [6}. 
However the latter is much smaller for all but the largest axigluon masses and we will ignore this process for the 
discussion in this Section. 
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FIG. 4: Axigluon pair production cross section from gluon initial states at the 7 TeV LHC as a function of 
the axigluon mass for x = 1 (solid) and x = (dashed). Also shown are the upper bounds on this cross 
section obtained by ATLAS (HIE! and CMS EH which apply if axigluons decay with 100% branching 
fraction to dijets. The latter has been unfolded by comparing the cross-section times acceptance for col- 
oron (axigluon) pair production presented in the analysis with our axigluon pair production cross-section 
calculated with MadGraph/MadEvent 5. 

analyses at ATLAS and CMS. The bounds apply only to the case where the axigluon decays 
predominantly to dijets. The two ATLAS analyses focus on events with 4 hard jets and form 
the invariant masses of all possible combinations of jet pairs. Keeping only events in which the 
invariant masses of the two pairs of dijets are similar to each other one can enhance signal over 
background sufficiently. The 2010 ATLAS [|29l analysis with 36 pb _1 takes advantage of the 
low instantaneous luminosity and low jet p T triggers of the 2010 LHC run to search for light 
dijet resonances in the mass range from 100 to 200 GeV. The recent analysis of 4.6 fb _1 of data 
from 2011 11301 uses a similar technique with higher jet thresholds. It is sensitive to axigluon 
masses from 150 GeV to 350 GeV. Finally, the CMS analysis OTTl looked for pair-produced dijet 
resonances in events with 4 hard jets (E T > 150 GeV) and employed a "diagonal cut" ll53l in 
order to enhance signal to background. The CMS analysis is sensitive to axigluon masses from 
320 GeV to 580 GeV. Each of the three analyses excluded cross sections well below the axigluon 
one, and combining limits from all three analyses covers the axigluon mass range from 100 GeV to 
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580 GeV. Thus the simplest version of the light axigluon model with no new light particles below 
the axigluon mass is ruled out even when allowing for flavor non-universal axigluon couplings to 
quarks. 

In the next Section we will discuss how these constraints can be evaded by opening up new 
axigluon decay channels with large branching fractions into multiple jets. This suppresses the 
branching fraction to dijets and creates new possibilities for discovering the axigluon in multijet 
resonance searches. 

For completeness we now list other constraints on axigluon models which are independent of 
the above mentioned dijet resonance pair searches. However we emphasize that the dijet resonance 
pair searches alone are sufficient to rule out axigluons decaying to dijets. 

• Dijets: Axigluons can be singly produced as an s-channel resonance decaying into quark- 
antiquark pairs. If the axigluon is sufficiently narrow it would lead to a clearly visible 
resonance in the dijet invariant mass spectrum. Assuming that there are no new decay chan- 
nels for the axigluon, and assuming flavor-universal couplings to quarks the axigluon width 
rU = 5/(24^)^771^ is negligible compared with experimental resolution. Then the di- 
jet searches by UA2 ll54ll and CDF [55J rule out axigluons with the necessary couplings 
to explain the tt asymmetry for all masses from 140 to 400 GeV. The dijet bounds can be 
evaded by reducing the couplings to first generation quarks and increasing the coupling to 
top quarks. We discuss dijet bounds in more detail in Appendix [Bj 

• Precision low energy measurements: Axigluon couplings to quarks are also constrained by 
loop corrections to the Z coupling to quarks. These constraints are only significant for the 
smallest axigluon masses near 100 GeV OTlBTi . Another potential constraint [|28l derives 
from axigluon loop corrections to the Z coupling to quarks as measured in atomic parity 
violation. The couplings in our model are too small to give a significant effect. 

• Decays to b-quarks: In the case of flavor non-universal axigluon couplings one might expect 
enhanced branching fractions to bottom quarks. CDF ||56l has performed a search for a 
similar final state motivated by Higgs production in association with additional b-quarks. 
Assuming 100% branching fraction of axigluons to b quark pairs we find that the CDF 
study can be used to rule out masses up to 250 GeV. 
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FIG. 5: Diagrams contributing to axigluon pair production from initial state gluons: (a) gluon s-channel, 
(b) axigluon t-channel, (c) 4-point interaction. 
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FIG. 6: Diagrams contributing to axigluon pair production from a qq initial state: (a) gluon s-channel, (b) 
SM quark t-channel, (c) heavy quark t-channel, and (d) axigluon s-channel. All four diagrams have to be 
included with appropriately chosen axigluon couplings to preserve unitarity. 

IV. NEW AXIGLUON DECAY CHANNELS AND MULTIJETS 



In this Section we explore the possibility that the axigluon is broad because it has new decay 
channels with large partial widths. The cross section for pair production of axigluons of mass 
m,A is only about an order of magnitude smaller than the huge QCD dijet cross section with jet 
transverse momenta pt > txia- Therefore these new decay channels have to be very difficult to 
detect or buried in QCD background to not have been ruled out already. This eliminates significant 
branching fractions to final states with leptons, photons or missing energy, and leaves decays to 
jets as the only realistic possibility. Figure|7]shows the cross section for axigluon pair production at 
the LHC with 8 TeV center of mass energy. It is dominated by the model independent gg — > AA 
process, with the qq initiated process contributing only about 10%. If the qq — > AA process is 
included in event generations care must be taken to employ a consistent unitary model including 
all diagrams in Figure|6](see Appendix|C]). For all axigluon masses between 100 and 400 GeV the 
cross section at the LHC is very large, producing 10 9 — 10 6 events in the 2012 run. Thus if the 
signal can be isolated from QCD backgrounds the axigluon should be observable in multijets. 

In the previous section we discussed axigluon decays to dijets and showed that axigluons which 
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FIG. 7: Axigluon pair production cross section at the 8 TeV LHC. Shown is the cross section with x = 1 
(solid) and x = (dashed). In both cases the cross section is dominated by the process gg — > AA but 
includes the smaller qq — > AA. The cross sections have been computed in the parity symmetric model of 
Appendix |A| with heavy quark masses set to 150 GeV, gA = and g m i xe d = g s - 

decay predominantly to dijets are ruled out over the entire mass range of interest. Here we focus on 
axigluons which decay preferentially to 3, 4, or 6 jets, giving events with 6-12 jets from axigluon 
pair production. Discovering or ruling out the axigluon is a matter of systematically eliminating 
the possible decay topologies. To this end we define four "simplified models" intended to study 
the four most "reasonable" decay topologies. The models are distinguished by different axigluon 
decay topologies which result from different intermediate particles through which the decay can 
proceed. We define the couplings of the new particles and the axigluon as coefficients in an effec- 
tive Lagrangian. Finally, we discuss collider signatures of each model and point out where existing 
searches already limit the allowed parameter space. For simplicity and for ease of comparison we 
assume 100% branching ratios to the selected final states. MadGraph implementations for each of 
the models are provided here 041 . UV-complete models which serve as explicit examples for the 
scenarios described here and demonstrate their consistency are presented in Appendix |A~} 

Finally, before we discuss the models and the signatures of pair production of axigluons at the 
LHC we comment on single production of axigluons at the Tevatron via qq — > A. If the axigluon 
predominantly decays to 3, 4 or 6 jets as postulated here one should also be able to see it as a 
resonance in 3, 4, or 6 jet events at the Tevatron. While a simple multi-jet resonance search is 
straightforward (one looks for a resonance in the total invariant mass spectrum of multijet events) 



13 



we are not aware of any public results 5 except for a preliminary analysis from DO 115711 using 
0.7 fb _1 and focusing on resonances above 400 GeV. For the case of flavor-universal axigluon 
couplings to quarks we expect that if any such resonance searches were performed at the Tevatron 
they could rule out axigluons in all of our models, at least for the heavier masses that we consider. 
For the lightest masses of order 100 GeV the jet-pr thresholds at the Tevatron may already be 
too high. We urge that such multijet resonance searches at the Tevatron will be carried out and 
published as they would lead to strong limits. However, by allowing non-universal couplings it is 
possible to reduce the coupling to light quarks and reduce the single axigluon cross section by as 
much as two orders of magnitude. This would allow the multijet resonance to hide in the QCD 
background. 

In the following, we return to discussing pair production at the LHC (and Tevatron). It is much 
more difficult to identify the resonances in pair production due to combinatoric backgrounds. 
However pair production has the big advantage that the cross section is model independent with 
no parameters which can be tuned to suppress it. 

Scenario A, eight jets with scalar intermediate resonances 



A" 




FIG. 8: Axigluon decay to two scalar (or pseudoscalar) color-octets with subsequent decay to four gluons. 

In this scenario there is a color-octet of scalars n = n a T a with masses < with strong 

couplings to axigluons so that axigluons decay predominantly to pairs of these scalars. The n a 

5 The dijet search from CDF |55| includes multijet events in the data because the analysis does not veto additional 
jets. However only the two leading jets are used to form the invariant mass which leads to significant smearing and 
reduction of the reconstructed masses, and no significant bounds can be obtained. 
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further decay to pairs of gluons. Thus axigluon pair production leads to final states with eight jets 
with the resonance structure shown in Figure [8] This model is nice in that it has a simple signature 
and it is easy to construct consistent models with the necessary couplings. A drawback is that the 
introduction of the color-octet scalars is ad hoc. The scalars of this model can also be searched for 
directly through QCD pair production with subsequent decay to 4 jets. Existing searches do not 
exclude masses below 100 GeV down to 10s of GeV if the scalar carries no electroweak quantum 
numbers and has only small couplings to quarks. Above 100 GeV only a small mass window 
around 140 GeV is allowed by the aforementioned ATLAS search OOl . 

Lagrangian 

The Lagrangian describing the couplings of the scalar to axigluons and gluons G a is 

c = \ {D^r (D^r + x A r bc A a y(D^r + -4- tr (7r G G ^ } (3) 

where (D^n) a = + g s f abc G b ^K c and = d^G v - d v G^ - ig s [G^, G„] is the gluon field 
strength. The coupling constant must be at least of order 1 so that axigluon decay to scalars 
dominates over decay to quarks. The mass parameter A determines the scalars' decay width to 
gluons. In a UV completion it may correspond to the mass of new vector- like fermions, and the 
operator tr (% G^G^) is generated when integrating out these fermions. 6 

Phenomenology 

How the multi-jet signal materializes depends on the relative mass of the axigluon and the 
scalar. In order for the axigluon decays to scalars to dominate over decay to quarks the scalar 
must be significantly lighter than half the axigluon mass. On the other hand, for m n ^ m^/8 the 
opening angle between the two gluon jets from each scalar decay becomes so small that the jets 
merge. Then the 8 gluon final state would be observed as 4 jets and the CMS and ATLAS bounds 
in Figure [4] apply. Thus for this model to be viable scalar masses must lie roughly between 1/8 
and 1/2 of tua- 

6 Several modifications of this simple model are possible. Instead of the scalar we could have a pseudoscalar which 
couples to gluons as tr ^7r G^G^J . The scalar or pseudoscalar may also decay to quark anti-quark pairs, and one 
would expect the heaviest accessible quarks, the b-quarks, to dominate. 
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Given the large cross section there are several options for probing this model. Perhaps the most 
straightforward is to look for a threshold feature at twice the axigluon mass in final states with 
large numbers of jets. An alternative would be to make use of the resonance structure of the events 
and identify resonances in 2 or 4 jet invariant mass distributions corresponding to the scalars and 
axigluons. 

A search of the former kind which could potentially be sensitive to axigluon pair production is 
the CMS Black hole search [58]. Here one looks for an enhanced rate of multi-jet final states at 
large invariant masses or large St- The publicly available analysis is not sensitive to axigluon pair 
production near threshold for the axigluon masses of interest because of high jet-pr thresholds 
and a focus on very large St > 1-8 TeV. However one might still see an enhancement of the 
cross section at the highest energies due to the large tail of the axigluon pair production cross 
section. This is more difficult than looking for the threshold bump because it requires a theoretical 
prediction of the QCD background. Another issue is that the number of jets observed depends on 
the relative size of jet-pr cuts, the axigluon and scalar masses and typically some jets are lost. On 
the other hand, the number of jets can also exceed 8 due to the large amount of QCD radiation 
expected from an 8 gluon final state. 

The alternative search strategy would be to look for invariant mass bumps in multijet events. In 
principle one can try to identify both the mass of the axigluon as a 4 jet resonance and the mass 
of the scalar as a 2 jet resonance. 7 The 2 jet resonance is easier to reconstruct, especially if the 
scalars are boosted so that the two gluon jets from their decay are near each other. For example, 
if rriA/4:, then one might define broad jets (R ~ 1) to capture the two gluon jets from 

each scalar decay in one broad jet. Then jet-substructure techniques ll60ll could be used to look 
for subjets within the broad jets which combine into a resonance at m w . The axigluon resonance 
is more difficult to reconstruct in general because of large combinatoric backgrounds. Axigluon 
bump hunting might be most promising in events where the axigluons are produced significantly 
above threshold so that the decay products of the two axigluons are boosted and do not overlap. 

For previous phenomenological work on a very closely related model with explicit proposals 
for cuts to enhance the resonance signals see for example ll52~ll6TTl . 



7 Simultaneous with the posting of this paper, CMS released a preliminary analysis searching for 8jet finals states 
with this topology which is sensitive to (and rules out) axigluon masses above 400 GeV (59]. Going to lower masses 
is more difficult due to trigger thresholds. 
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Scenario B, 12 jet final state 




FIG. 9: Axigluon decay to two heavy color-triplet fermions with subsequent decay to six quarks. 

In this scenario we consider new fermions (heavy quarks) which transform as color triplets like 
the ordinary quarks. Axigluons are assumed to predominantly decay into pairs of these heavy 
quarks. The heavy quarks subsequently each decay to 3 light quarks via off-shell axigluons. Thus 
in this model pair production of axigluons leads to 12 jet final states as shown in Figure [9} This 
model is natural in the sense that heavy vector-like quarks are already required by unitarity in a 
renormalizable axigluon model (see Appendix [A]). Unitarity and bounds from LEP2 constrain the 
masses of these heavy quarks to lie between 100 GeV and approximately 1 TeV. 8 For definiteness, 
we assume that only one of the heavy quarks is lighter than 1/2 of the axigluon mass so that only 
this fermion is involved in axigluon decays. All other fermions are assumed to be near 1 TeV 
so that they are irrelevant to axigluon decays. Alternatively, one could have also chosen multiple 
heavy quarks to be light. 

In principle there are constraints from direct pair production of the heavy quarks at the LHC 
with subsequent decay to six jet final states. However, the pair production cross section of the 
quarks is sufficiently small that even several flavors of such color-triplet fermions decaying to 
three jets are allowed by all searches 116211631 . 

8 We discuss the relationship between heavy fermion masses and unitarity of the axigluon pair production cross 
section in Appendix |C| 
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Lagrangian 

Using two-component spinors for the fermions, the Lagrangian describing the couplings of the 
axigluon to the vector-like heavy quark (D H , D H ) is 

C = X H A;D^T a D H + \ H A«D^T a *D H + (\ mix A*D^T a D SM + h.c.) . (4) 

Here D SM is a right-handed down-type SM quark, D H has the same quantum numbers as D$m 
and has a mass with its vector partner, M H q D h D h . Since we are only including a single heavy 
quark, the mixed coupling to Dsm necessarily breaks the SM flavor symmetries. Flavor constraints 
on the couplings of 1st and 2nd generation quarks are generally stronger than constraints on 3rd 
generation quarks. Therefore we chose the mixed coupling to only involve the right-handed bottom 
quark, i.e. D SM = b R . Then axigluon decays give final states with at least two b-jets. In order to 
ensure that the axigluon predominantly decays to pairs of heavy quarks one of the couplings \h 
or \h has to be larger than g s and also \ m i X -C max[A#, A#]. Explicit expressions for A#, Xh, 
and \ m ix in terms of parameters of a UV-complete model are given in Appendix [A] (and utilized 
in our MadGraph implementation of the model). These expressions imply relations between the 
couplings which ensure unitary amplitudes. 

Phenomenology 

Axigluon pair production in this model leads to a final state with even larger jet multiplicity 
than the one described in Scenario A. Combinatoric backgrounds to reconstruction of any of the 
resonances are therefore much larger. Thus any searches which focus on reconstructing axigluons 
or heavy fermions must rely on boosted events and possibly use jet substructure techniques. 

On the other hand, since the axigluon pair-production cross section is much larger than the 
background QCD 12-jet cross section one might expect to see a bump in the measured cross 
section for very high multiplicity jets. A search of this kind would be similar in spirit to the CMS 
black hole search [|58l but would have to use much lower p T cuts. There is significant sensitivity 
to jet-pT thresholds and the relative size of the axigluon and heavy quark masses, and in practice 
many of the jets may be lost due to cuts and jet merging. 
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FIG. 10: Axigluon decay to a heavy color-triplet fermion in association with a light quark giving rise to a 4 
quark final state. 

Scenario C, 8 jets with fermionic intermediate resonances 

In this scenario the axigluon is assumed to decay to a SM quark in association with a heavy 
vector-like quark. This heavy quark further decays to three light quarks so that axigluon decays 
lead to 4 light quark jets in the final state as shown in Figure [TOj In order to get a sufficiently large 
decay width to this final state we choose several such new heavy quarks. This scenario occurs 
naturally in renormalizable axigluon models in which the coupling of the axigluon to quarks is 
suppressed due to mixing with heavy quarks. In those models small axigluon couplings to SM 
quarks are correlated with large mixed couplings \ m ix of the axigluon to one SM quark and one 
heavy quark (see Appendix [A]). 

Lagrangian 

Using two-component spinors for the fermions, the Lagrangian describing the mixed couplings 
of the axigluon to the vector-like heavy quarks and light quarks is 

£ = A^ ix A*U H a»T a U SM + A£ ix A;D^T a D SM + h.c. . (5) 

Here U sm an d D SM are right-handed up- and down-type SM quarks, Uh (D h ) has the same quan- 
tum numbers as Usm (Dsm)- The heavy quarks have masses with vector-like partners Uh, Dh, 
for example MhqUhU h- Explicit expressions for A^{^ in terms of parameters of a UV-complete 
model are given in Appendix [A] (and utilized in our MadGraph implementation of the model). 
These expressions imply relations between g s , g A and A^f which ensure unitary amplitudes. 
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A few comments about the flavor structure of this scenario are in order. First, experimental 
searches for new heavy quarks with significant branching fractions to final states with W's and 
Z's rule out such new quarks with masses below 500 GeV. Mixing with third generation quarks 
generally introduces decays involving W's and Z's, and therefore we consider new quarks which 
only couple to 1st and 2nd generation quarks. This necessarily implies some level of flavor vio- 
lation. Constraints from flavor physics further forbid large couplings to quark doublets. This is 
why we did not include couplings to quark doublets in the Lagrangian in Eq.[5} Introducing heavy 
quark singlets without corresponding heavy quark doublets breaks parity symmetry in the axi- 
gluon coupling to heavy fermions. In a complete model one would therefore expect the axigluon 
couplings to light quarks to not be purely axial without fine tuning. 

To summarize, we consider two generations of degenerate heavy quarks, two of up-type Uh and 
two of down-type Dh and their vector-like partners. The couplings A^f are chosen to preserve 
the flavor symmetries of the 1st and 2nd generation. 

Phenomenology 

Axigluon pair production in this model leads to a final state with eight light quark jets as in 
model A. Most of the discussion for model A also applies in this case. The main difference is that 
because of combinatoric backgrounds reconstruction of a three jet resonance is more difficult than 
the reconstruction of two jet resonances. The other difference is that in this model there can be a 
significant asymmetry between the typical pr of the jets from the heavy fermion decays and the 
jet from the primary axigluon decay. This asymmetry depends on the relative mass of the axigluon 
and heavy fermions and contributes to the efficiency with which these events are picked up in 
searches similar to the black hole searches. 

Scenario D, interpolating between 6 and 8 jet final states 

This scenario is identical to scenario C except that there is an additional light color-singlet 
scalar particle r\ which participates in the heavy quark decay. Thus axigluons now decay through a 
small decay chain as shown in Figure [TT] Light singlet scalars arise naturally in axigluon models 
as uneaten components of the field(s) that break the UV gauge symmetry down to QCD giving 
mass to the axigluon. In the simple case where S77(3)i x SU (3)2 is broken to the diagonal 8 
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FIG. 1 1 : Axigluon decay to heavy color-triplet fermion in association with a light quark. The heavy quark 
subsequently decays to a SM quark and a light scalar which further decays to b-quarks. 

Nambu-Goldstone bosons are eaten to give a mass to the axigluon. However there is a 9th NGB, 
"axion", because the full global symmetry breaking structure is really Z7 (3)i x U(3) 2 broken to 
diagonal £7(3). This uneaten NGB is naturally light and can play the role of 77. Its mass can be 
chosen arbitrarily by adding small explicit symmetry breaking interactions. We envision 77 masses 
in the range from 10 GeV to the axigluon mass for this scenario. In the axigluon models described 
in |[T6ll the 77 is expected to have large mixed couplings with a SM quark and a heavy quark, so that 
heavy quarks preferentially decay into light quarks and 77. In addition 77 also has a small coupling 
to two SM quarks which arises from mixing proportional to the SM quark masses. This coupling is 
largest for the third generation quarks and causes the scalar to decay almost exclusively to bottom 
quark pairs. 

Lagrangian 

The Lagrangian describing this model is the same as model C, plus a new piece describing the 
axion interactions, 

C D ^(<9t7) 2 - ^77 2 + i\ b r](b[b R - b ] B b L ) + (\ a r)U H UsM + KvDhDsm + h.c.) . (6) 

In the equation above, [i a is the axion's mass, A a is the coupling that controls the decay of the 
heavy quark to a quark and the scalar and A& the coupling that controls the decay of the axion 
to the SM bottom quarks. \ b is not the SM bottom Yukawa coupling, but it is proportional to it 
because it is generated by mixing proportional to quark masses. 
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Phenomenology 

In this scenario the axigluon decays to a heavy quark and a light quark as in model C. The heavy 
quark then decays to a quark and the axion, which in turn decays to two bottom quarks. Thus each 



axigluon decays to 2 light quarks and 2 b quarks as shown in Figure 1 1 The phenomenology of 
this model is interesting as it interpolates between final states with 6 to 8 jets. When the rj mass is 
a sizeable fraction of the heavy quark mass all 8 jets are in principle observable as separate jets. 
However when the mass of r\ is closer to 10 GeV the 6-quarks from 77-decay become collimated 
and merge into a double-b jet. This double-b jet is quite interesting. It contains substructure due 
to the presence of two b-quarks and it is very likely to be tagged as a b-jet because it contains 2 
displaced vertices. Double-b jets are not unreasonable to occur in many models beyond the SM 
which contain light scalars and are an interesting signature to look for independent of the multi-jet 
final states we propose here. 

How would one look for this model? In the case of large 77 masses, the phenomenology is very 
similar to the 8 jet finals states which we discussed before. One can focus on searching for bumps 
in di-dijet invariant masses or take advantage of the high jet multiplicity of the events as discussed 
in Scenario A. For the smallest 77 masses the model is already mostly ruled out by the 6 jet final 
state searches for R-parity violating gluino decays ll62l l63l . If one takes the limit from these 
searches at face value there is only a small window in axigluon masses between 170-200 GeV 
(near the top mass) where this scenario is not ruled out. However the acceptance of the search 
to the axigluon decay signal depends on the mass of 77. Thus this model motivates an interesting 
possible extension to the current 6-jet resonance searches for R-parity violating gluino decays. 
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Appendix A: UV-complete axigluon models 

In this Appendix we present two complete models for the axigluon and discuss their connection 



to the simplified models described in section IV In both models the axigluon arises from the 



breaking of a larger gauge group 577(3) x 577(3) to the diagonal 577(3), which corresponds to 
the QCD gauge group. The first model has a parity symmetry built in under which the axigluon 
is odd. Parity ensures that the light quarks have axial couplings to the axigluon. In this model all 
particles couple to the axigluon with couplings that are bounded by g s . In the second model there 
is no parity symmetry and getting axial couplings to the axigluon requires fine tuning. On the 
other hand it is easy to introduce particles with large couplings to the axigluon. This is desirable 
because if these particles are light then the axigluon can have a large partial width to decay to them 



as required in models A and B in Section IV 



1. Parity symmetric two site model 



Q,U,D 
U',D' 





FIG. 12: Moose diagram for the parity symmetric two site model. 



This model can be described by the diagram in Figure 12 There are two SU(3) groups with 



equal gauge couplings. There is a scalar field <p which is a fundamental under 5£7(3) l and an anti- 
fundamental of SU(3)r. The Q is a left-handed Weyl fermion transforming as a fundamental of 
SU(3)l, and it has the same electroweak quantum numbers as the SM quark doublets. The fields 
U and D are right-handed Weyl fermions, transform as fundamentals of SU(3)r and have the 
electroweak quantum numbers of the up-type and down-type SM quark singlets. In addition, there 
are vector-like partners for each of these fields, e.g. (U', U) are the partners of U . The primed 
fields are fundamentals under the opposite SU (3) and have the same chirality and electroweak 
quantum numbers as their unprimed partners. The barred fields have the same chirality and the 
opposite gauge quantum numbers as the primed fields so that, for example, U' and U can have a 
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mass term MUU' . 

This model has a parity symmetry which corresponds to flipping the diagram in Figure[l2j e.g., 
Q <H- (U,D) and Al -h- Ar. Clearly this cannot be an exact symmetry because Q and (U,D) 
have different electroweak quantum numbers (the same reason why parity is not a good symmetry 
of the SM). Nonetheless, we assume that it is a symmetry of the extended strong interactions and 
corrections due to the weak interactions are small (suppressed by a loop factor). 

The gauge groups are broken to the diagonal by the VEV of the scalar, (0) = f Izxs- This 
gives a mass tua = \/2gf to the anti-symmetric combination A^ = i (A^ — A R ) and leaves the 
symmetric combination massless. The massive vector A^ is the axigluon and the massless one the 
gluon G^. Fields charged under SU(3)l/r couple to the axigluon with couplings ±g s = ±g/\/2. 

In order to get suppressed couplings of the SM quarks to axigluons we now introduce mixing 
between the fermions. Since the fermions and their partners have opposite sign couplings to the 
axigluon any linear combination of the two will have a reduced coupling. The mixing is obtained 
by adding a Yukawa coupling involving the link field <p in addition to the Dirac mass of U' with 
U, e.g. 

C = U (MU' + \<f>U) + h.c. (Al) 
After plugging in the VEV for <f> one sees that the combination Uh = cos aU' + small gets 



a mass Mhq = a/ M 2 + A 2 / 2 with U, where tana = Xf/M. The orthogonal combination 
Usm = cos aU — sin all' remains massless and corresponds to the SM quark; it eventually gets a 
mass from the SM Higgs VEV via Yukawa couplings which we have not explicitly displayed here. 
The couplings of the mass eigenstates to the axigluon are given by 



C = g s A a cos(2a)U j H a f "T a U H - cos(2a)^ M T a f/ SM - sm{2a){U ] H a^U SM + h.c.) 



. (A2) 



We see that the axigluon coupling to SM quarks is g& = g s cos 2a. By choosing appropriate values 
for a we can obtain — 0.3 — 0.5 as required to explain the tt asymmetry. The couplings of the 



axigluon to SM quarks are automatically axial as long as the mass and Yukawa terms Al for the 
left-handed, Q and right-handed fermions (U, D) and their respective partners respect the parity 
symmetry. The axigluon also has mixed coupling allowing transition between a SM quark and its 
heavy quark partner with g mix = — sin(2 a)g s , the couplings satisfy the relation g 2 s = g\ + g^ix- 
To summarize, this extension of the SM contains an axigluon with mass 2g s f. The SM 
quarks get their masses from the Higgs VEV as usual and couple axially to the axigluon with 
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coupling gA = g s cos(2a). The model contains heavy partner quarks 9 with masses Mhq = 
a/M 2 + A 2 / 2 . The axigluon has mixed couplings to SM quarks and their heavy partners with 
9mix = —9s sin(2a). It also couples to two heavy quarks with couplings —g& (to U H ) and g s (to 
U). 



Connection with the simplified models of Section IV 



The model discussed in this section provides a good skeleton for building complete versions of 



the simplified models C and D presented in Section IV The coupling A m between the axigluon, a 
quark and its heavy partner is fixed to be g mix = —g s sin(2a). Because this coupling is bounded by 
g s one needs multiple flavors of heavy quarks being lighter than the axigluon in order to generate 
a large enough partial decay width of axigluons to quarks and their heavy partners. 

As discussed in the Lagrangian description of scenario C we choose all partners of the 1st and 
2nd generation right-handed quarks to be degenerate and light, this corresponds to 4 heavy quarks 
lighter than the axigluon. Notice that this explicitly breaks the parity symmetry in the axigluon 
sector, since we are treating the partners of U's and D's differently than the partners of Q's. This 
leads to some fine-tuning in order to preserve the purely axial couplings of the axigluon to SM 
quarks. Models with more sophisticated flavor structure than the one presented here might allow 
one to preserve the parity symmetry. We will not pursue this issue further in this paper. 

If the "axion", i], which is a component of the link field (ft, is given a mass larger than that of 
the lightest heavy quark than one has an implementation of model C. If the axion is lighter than 
the heavy quark one has an implementation of model D, with the couplings of t] to SM quarks 



induced by the Yukawa interaction in Eq. Al 



2. Asymmetric two site model 



This model can be described by the diagram in Figure 13 The axigluon also arises from the 



breaking S77(3)i x S77(3) 2 — > SU(3) Co \or- However in this case the gauge couplings are not 

9 As it stands this model has gauge anomalies. The anomalies can easily be cancelled with additional chiral fermions 
which are vector-like under the SM gauge group. Their mass must be proportional to the <fi VEV and is bounded 
by 47r/. Therefore these fermions can be pair-produced at the LHC. For simplicity, we assumed that they are too 
heavy to play a role in axigluon phenomenology. 




FIG. 13: Moose diagram for the asymmetric two site model. 

equal and thus the linear combinations corresponding to the gluon and the axigluon are no longer 
symmetric and anti- symmetric in A\ and A 2 . Defining tan f3 = g/G we have 

G" = cos /3A% + sin f3A% , A^ = sin /3A% - cos f3A% . (A3) 

The axigluon mass is itla = a/(? 2 + G 2 f, and the strong coupling is given by g s = gGj \/ 'g 2 + G 2 . 

In this model there are left-handed fermions Q and right-handed fermions U, D all charged 
under S77(3)i. They end up with couplings to the axigluon given by ^g s which can be made 
smaller than g s simply by taking g < G. This is nice because small axigluon couplings are 
needed to explain the ti asymmetry. The problem is that the couplings are purely vectorial. This 
can be fixed by introducing fermion mixing. For definiteness we choose to only mix the right- 
handed fermions U, D with right-handed heavy fermions U', D' charged under SU(3) 2 . The 
massless combinations which are identified with the SM fields are Usm = cos all — sin all' and 
Dsm = cos aD — sin aD'. Their couplings to the axigluon are given by 

9a = 9s(^ cos 2 a - ^ sin 2 a^j , (A4) 

where a is the mixing angle between unprimed and primed fields. We see that by fine-tuning the 
mixing angle so that sin a = V%g/ a/ g 2 + G 2 = \/2g s / G one can obtain axial couplings between 
the axigluon and the SM quarks. 

Notice that any particle which is charged only under SU (3) 2 couples to the axigluon with cou- 
pling — g s G/g. Therefore one can easily introduce particles that couple strongly to the axigluon 
by making them charged under SU (3)2. This allows for a large decay width of the axigluon to 
such particles if they are lighter than the axigluon. 
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The Lagrangian describing the couplings of the axigluon to the gluon is then 

C = -hi(F^F, u ) + m 2 A tv(A^) + ig s tr (G^A", A"]) (A5) 

+ Ws (§ - |) tr {F^W, A«\) + f (l + (G?2 g ~f )2 ) tr A v f) , 

where is the gluon field strength and F^ u = D^A„ — D V A^ is the axigluon field strength. 
The third term in the equation above is the x-term of Eq. [2] with x = 1 as required by unitarity. 
The fourth term contains a triple axigluon vertex and is absent in the parity symmetric limit when 
g = G. 

This model is another example of a UV completion for the axigluon model. The mass of the 
axigluon is given by niA = \/ g 2 + G 2 f. In addition to the SM quarks there are heavy vector-like 
partners for all right-handed SM quarks (note that this was an arbitrary choice, there could instead 
be partners for the left-handed quarks or for both). The couplings of SM quarks to the axigluon 
are naturally suppressed, but in order to obtain axial couplings one needs to fine tune the quark 
mixing angle a. Despite this ugly fine tuning this model has a few advantages over the previous 
one: it requires a smaller number of extra particles, one can easily introduce decay channels with 
large partial widths for the axigluon, and the fermion assignments are automatically anomaly free. 



Connection with the simplified models of Section IV 



Using the construction described in this section one can easily implement a complete version 



of models A and B of Section IV which require large couplings of new intermediate resonances 
to the axigluon. To implement Model B we take the mass of the heavy partner of the bottom quark 
to be smaller than half the axigluon mass and all other heavy quarks heavier than the axigluon. We 
also choose sin a = y/2g s /G The couplings of the axigluon to the heavy quarks defined in Eq.|4] 
are then 



— 9s— [~Q2 ~ ) ' Xh = 9sj, Knx = g s y 2 - — . (A6) 



We see that one automatically has a large coupling (since we are taking G > g ) of the axigluon to 
two heavy quarks as required in model B. 

In order to implement model A one assumes that all heavy quarks are heavier than the axigluon, 
so that it cannot decay to heavy quarks. Then one includes an additional scalar n which is an ad- 
joint of SU (3) 2 and not charged under SU (3) i. After the breaking SU(3) l xSU{3) 2 -»■ SU(3) Co i t 
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the scalar becomes an adjoint of the QCD gauge group as described in model A. Because it is 
charged under SU(3) 2 one finds that it's coupling to the axigluon, A a of Eq.[3j is given by —g s G/g 
and thus is automatically enhanced compared to the QCD coupling. In order for the scalar to decay 
we introduce the dimension 5 operator tr {txG^ u G^ v ) which couples the scalar to two gluons. This 
operator may be generated by integrating out an additional vector-like pair of fermions charged 
under SU (3)2 which have Yukawa coupling to n. These fields could be very heavy so that their 
only phenomenological consequence is the dimension 5 scalar decay operator. 



Appendix B: Constraints from axigluons as dijet resonances 

In this Appendix we consider dijet constraints 10 on axigluon models from UA2 and Tevatron. 
Since the axigluon can be produced from a qq initial state it can also decay into qq giving rise to 
dijet events. Assuming that the axigluon couplings to quarks are flavor-universal (g A = g A ), the 
dijet cross section depends on the same parameters tua, ou, as the tt asymmetry and strong 
constraints can be obtained. In particular, from Table [I] and Figure 14 we see that couplings g& in 



the range 0.3 to 0.45 lead to the desired 10% tt asymmetry from NP. For such small couplings the 
axigluon width Ya = B/(24:ir)g A mA is always negligible (the experimental resolution is on the 
order of 10%) and narrow resonance searches apply. In Figure [14] we compare the limits obtained 
from dijet resonance searches by UA2 [|54l and CDF [1551 to contours of constant tt asymmetry. 
To obtain these limits we computed the axigluon-mediated dijet cross sections to leading order 
using MadGraph and compared with the cross section limits quoted by the experiments. We see 
that over the whole range of masses where searches are available (140-400 GeV) this simplest 
axigluon explanation of the tt asymmetry is inconsistent with dijet constraints. One also sees that 
the limits are relatively weak for axigluon masses near 280 GeV. Until recently this would have 
motivated flavor non-universal axigluon models with reduced couplings to first generation quarks 
and enhanced couplings to top quark to compensate. However such models have been ruled out 
by LHC searches for pair production of heavy resonances which decay to dijets as discussed in 
Section Iml 

10 For a recent comprehensive review of dijet constraints from hadron colliders see 
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FIG. 14: Contours of constant tt asymmetry from axigluon exchange (thin black) versus upper limits on the 
axigluon coupling gA from dijet searches at UA2 (dashed red) and CDF (solid red). 

Appendix C: Unitarity 



In this Appendix we discuss unitarity constraints on axigluon couplings and demonstrate po- 



tential pitfalls with an explicit example (Fig. 15 ). Axigluons are massive vector bosons, and in a 
weakly coupled theory they must arise from a spontaneously broken gauge symmetry. The broken 
symmetry imposes relations between the coefficients of different terms in the Lagrangian which 
ensure cancellations between different diagrams when computing scattering amplitudes. These 
cancellations are required to prevent scattering amplitudes from becoming unphysically large at 
high energies and spoil unitarity. This is very familiar from the SM where the e + e~ — > W + W~ 
scattering amplitudes from individual diagrams diverge at high energies but are well-behaved once 
summed together. 

In multijet searches unitarity constraints can become especially relevant because many analyses 
require hard cuts to jet energies in order to suppress QCD backgrounds. These analyses are only 
sensitive to the high energy tail of axigluon pair production. If one uses an inconsistent model in 
which unitarity is violated or leaves out some diagrams in the computation of the axigluon signal 
this tail can be overestimated by orders of magnitude. 

At the LHC, to leading order there are two independent axigluon pair production modes. First, 
axigluons can be produced from a gg initial state as in Fig. |5J As discussed in Section III. this 
cross section is unitary as long as one includes the contribution from the tr (G^ U [A^, A,,]) -term 
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FIG. 15: Leading order differential cross sections for pair production of a 200 GeV axigluon as a function 
of axigluon px at the 7 TeV LHC. The black curve corresponds to a consistent model that includes 150 GeV 
heavy quarks partners of 1st and 2nd generation quarks and with X = 1. The red dashed curve corresponds 
to a model that includes the 150 GeV heavy quarks but with \ = 0- The blue dotted curve corresponds to 
a model without heavy quarks and with % = 1. In all 3 models we set qa = and the coupling between a 
quark and its heavy partner to be equal to g s . 

in Eq. [2] with x — 1- This term arises automatically in UV-completions of the theory, e.g. in the 
models of Appendix |A| and the two and three site models discussed in lfT6ll . 

The second axigluon pair production mode at the LHC is from quark-antiquark collisions as in 
Fig. |6j At the 7 or 8 TeV LHC this mode is much smaller than the gg initiated mode except for 
events with very high invariant mass, when the qq parton luminosities become larger than the gg 
luminosities. Nonetheless the qq initial state can be very important because many experimental 
analyses impose hard cuts suppress QCD background. Phenomenological studies for such analyses 
must employ simulations with both gg and qq initial states, and one must be careful to use a model 
with consistent quark couplings to axigluons. As we explain in the following paragraphs this 
requires the existence of new heavy fermions, and their contributions to axigluon pair production 
at high energies cannot be ignored. 

We showed in section III that the coupling of axigluons to SM quarks must be less than the 
QCD gauge coupling g s in order to explain the ti asymmetry. Small axial couplings of axigluons 
to quarks can be obtained through fermion mixing thus requiring heavy fermion partners (to mix 
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with). These heavy fermions contribute to axigluon pair production through t-channel exchange 
as shown in Fig. [6] (c). Neglecting these contributions produces unphysical cross-sections which 
violate unitarity. Moreover, in the absence of a parity symmetry there is also a contribution from 
a vertex with 3 axigluons as shown in Fig.[6](d). 

In Fig. [15] we show the cross sections for axigluon pair production at the LHC for the parity 
symmetric model of Appendix [A] (black curve) and for two inconsistent models (blue and red 
curves). The plot illustrates the unphysical growth of amplitudes in inconsistent models. The blue 
curve corresponds to a model with x = 1 but without heavy partners for the light quarks. Hence 
the contribution from the diagram in Fig. [6](c) is missing, and there is a large enhancement of the 
partonic cross section qq — » AA at high energies (the cross section falls after convolution with the 



parton distribution functions but the amplitude actually grows in this model). In Fig. 15 this results 
in an increase by more than a order of magnitude of the axigluon pair production cross section at 
high p T . The red curve corresponds to a model that includes heavy quark partners consistently but 
which has x = 0. In this model both the gg — > A A and qq — > A A amplitudes violate unitarity 
at high energies. In addition, the coupling proportional to x a ls° has an important effect on the 
gg — > AA amplitude at low energies so that setting x = results in a significant underestimate of 
the cross section near threshold. 
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